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Abstract. Cooling of neutron stars with dipole magnetic 
fields is simulated using a realistic model of the anisotropic 
surface temperature distribution produced by magnetic 
fields. Suppression of the electron thermal conductivity of 
outer stellar layers across the field increases thermal isola- 
tion of these layers near the magnetic equator. Enhance- 
ment of the radiative and longitudinal electron thermal 
conductivities in quantizing magnetic fields reduces ther- 
mal isolation near the magnetic poles. The equatorial in- 
crease of the isolation is pronounced for B > 10 10 G, while 
the polar decrease - for B > 10 12 G. The effects compen- 
sate partly each other, and the actual influence of the 
magnetic fields on the cooling is weaker than predicted by 
the traditional theories where the equatorial effects have 
been neglected. 

Key words: stars: neutron - dense matter - magnetic 
fields - radiation mechanisms: thermal 



1. Introduction 

Thermal X-ray radiation has been observed recently with 
ROSAT from at least four cooling neutron stars: PSR 
0833-45, PSR 1055-52, PSR 0656+14 and Geminga (see, 
e.g., Ogelman 1995). A soft (T s - (5 - 10) x 10 5 K) com- 
ponent of the observed spectra is most probably emit- 
ted from all the visible neutron star (NS) surface. It 
pulses with the pulsar period, and the pulsed fraction, 
~ 10 — 30%, varies with photon energy. The spectra and 
the light curves of these middle-age (10 4 - 10 6 yr) NSs 
contain important information on NS parameters and on 
the properties of superdense matter in NS interiors. 

The pulsations are most probably caused by the tem- 
perature variation over the NS surface. The natural cause 
of the variation is non-uniformity of the surface magnetic 
field and anisotropy of the heat transport in strongly 
magnetized subphotospheric layers of cooling NSs (e.g., 
Yakovlev & Kaminker 1994). 

So far most of the cooling calculations of NSs (e.g., 
Nomoto & Tsuruta 1987, Van Riper 1991) have been per- 



formed under the traditional simplified assumption that 
the magnetic field is radial everywhere over the stellar 
surface. In this case the thermal energy is carried through 
the subphotospheric layers along the magnetic field lines, 
and the surface temperature T s is uniform and determined 
by the longitudinal (along the field B) thermal conduc- 
tivity. The latter conductivity is mainly enhanced by the 
magnetic field, which increases T s at given internal tem- 
perature T\. As a result, the traditional cooling theories 
predict (e.g., Van Riper 1991) that the strong magnetic 
fields B > 10 12 G increase T s (and the photon luminos- 
ity) at the neutrino cooling stage. At this stage, a NS is 
not too old (age t < 10 4 -10 5 yrs) and cools mainly via 
neutrino emission from the stellar interior. On the other 
hand, strong fields decrease T s and accelerate the cooling 
at the subsequent photon cooling stage when the neutrino 
emission becomes low and the star cools via the thermal 
surface radiation. 

The main disadvantage of the traditional approach is 
that it neglects those parts of the NS surface where B 
is essentially non-radial, and the transverse thermal con- 
ductivity is important. The aim of this paper is to use a 
realistic model of heat transport and the related surface 
temperature distribution valid for any magnetic field ge- 
ometry in the NS surface layers. The model generalizes the 
results of Greenstein & Hartke (1983), Van Riper (1989), 
Schaaf (1990a), and Page (1995) who performed detailed 
studies of the NS surface temperature either for restricted 
magnetic field geometries or under specified assumptions 
on thermal conductivity of stellar matter (see Sect. 2, for 
details). Adopting this model we will carry out the cool- 
ing calculations (Sect. 3). Our results show (Sects. 4 and 
5) that the effects of the magnetic fields on the NS cooling 
are more sophisticated than anticipated previously. 

2. Surface temperature distribution 

Consider the temperature distribution over the surface of 
a magnetized NS. Let the star be not too young (t > 
1 - 10 3 yrs), so that the internal thermal relaxation is 
achieved. Then the stellar interior is highly isothermal: 
T = T(r) exp($(r)) is constant within the star, where 
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T(r) is the local internal temperature, r is radial coor- 
dinate, <J>(r) is the gravitational redshift function which 
takes into account the effects of General Relativity (e.g., 
Glen & Sutherland 1980). The isothermal interior is sur- 
rounded by a very thin subphotospheric layer. Its thick- 
ness ( ^ several meters) is much smaller than the stellar 
radius R ~ 10 km. This layer produces the thermal iso- 
lation of the NS interior. The temperature at the bottom 
of this layer is 7] = Texp(— <£> s ), where <fr s ss $(-R) = 
0.51n(l - R g /R), R g = 2GM/c 2 is the gravitational stel- 
lar radius, and M is the stellar mass. While the star 
cools down, the isolating layer becomes thinner (e.g., Gud- 
mundsson et al. 1983, Yakovlev & Kaminker 1994). The 
layer consists of two sublayers. In the inner one, the elec- 
trons are degenerate, and they are the main heat carriers. 
In the outer sublayer, the electrons are non-degenerate, 
and heat is mostly transported by radiation. In a rather 
cold NS, the surface layer may become solid, and the non- 
degenerate sublayer disappears. 

The effective surface temperature T s is determined by 
the heat transport through the isolating magnetized layer. 
The problem of calculating T s is complicated (Yakovlev 
& Kaminker 1994), and the exact solution has not yet 
been found. Let us discuss a realistic model solution. 
Since the isolating layer is thin it is reasonable to adopt a 
plane-parallel one-dimensional approximation almost ev- 
erywhere over the surface. We assume that the magnetic 
field is constant in each local part of the surface layer. 
Let the thermal flux F be normal (radial) to the surface 
and the tangential flux be insig nificant. Then F = crT s 4 , 
where T s is the local effective surface temperature, and a 
is the Stefan - Boltzmann constant. It is well known that 
T s depends on the internal temperature 7}, and on the 
surface gravity g s = {GM/R?)/^f{l - Rg/R) (T„ cx 5s 1/4 ; 
see, e.g., Gudmundsson et al. 1983). In our case, the local 
surface temperature T s — T s (B,x) depends also on the 
local magnetic field B and on the angle x between B and 
the normal to the surface. If T s (B,x) is known, the NS 
total thermal luminosity is given by 



L(T i )=aR 2 J T s 4 (B lX )d{l 



(1) 



(df2 being solid angle of a surface element), and the red- 
shifted luminosity (for distant observers) is L°° = L(l — 

Rg/R). 

Let the thermal energy in the isolating layer be carried 
by the thermal conduction (no convection or turbulent 
motion). Then the thermal balance yields 



(km cos 2 x + K -L sin 2 x) 



dT 



F. 



(2) 



where km and k± are, respectively, the longitudinal and 
transverse (with respect to B) thermal conductivities, and 
z is depth from the surface within the star in the lo- 
cally flat reference frame. The equation of thermal balance 



should be supplemented by the equation of hydrostatic 
equilibrium; we assume (as in all cooling theories) that 
the magnetic field in the NS surface layers is force-free. 
Both equations should be integrated from the surface in- 
side the star yielding the temperature and density profiles 
in the outer NS layer with given flux F (i.e., with given 
T s ). With increasing z, the temperature T(z) tends to a 
constant value to be defined as the internal temperature 
Tj. Combining the solutions for different T s , B and x one 
can determine the dependence of the effective local surface 
temperature T s on 7j, B and x- 

Since the thermal conductivities km and kj_ depend 
generally on T, B and density p the solution is complicated 
even with the above simplifications. Instead, we shall use 
the following model solution (Greenstein & Hartke 1983) 



T S \B, X ) = THB) cos 2 x + T±(B) sin 2 X , 



(3) 



where Tn(B) = T s (B,0) is the surface temperature for 
the case when B is normal to the surface (at the mag- 
netic pole), and T±(B) = T s (B,ir/2) is the same for the 
case when B is tangential to the surface (at the magnetic 
equator) . 

Let us discuss briefly the properties of Eq. (||). First, 
(||) yields exact solution of the thermal conduction prob- 
lem (||) for B normal or tangential to the NS surface. 
Second, (||) reproduces the limit of B = 0. Third, one 
can easily prove that (||) provides the solution for the case 
when Km k>±_ . The case is realized for not very weak 
magnetic fields and not too close to the magnetic equa- 
tor (k||COS 2 x 3> n±sm 2 x)- This case is very important 
since the transverse electron conductivity of degenerate 
electrons is much smaller than the longitudinal one, typ- 
ically, at B > 10 10 G (e.g., Yakovlev & Kaminker 1994). 
Fourth, Eq. (y|) is exact for the case when k_l/km is con- 
stant (independent of depth z) for a given local part of 
the NS surface layer. The latter case is also important. 
For instance, it is realized in a strongly degenerate non- 
relativistic electron gas at not too high fields B if the 
electron conductivity is produced by Coulomb scattering 
on atomic nuclei and the nuclei constitute a gas or liquid 
(e.g., Yakovlev & Kaminker 1994). The model solution (||) 
can actually be inaccurate at high magnetic fields near the 
magnetic equator where the thermal flux tangential to the 
surface can be substantial and the plane-parallel approx- 
imation can fail. However these areas of the NS surface 
give small contribution into the total stellar luminosity 
([!]), and, thus, they seem to be insignificant for cooling 
simulations. 

Using Eq. (||) one can obtain the following inequality 
of the thermal conduction problem 



T s (B,x) > T||(i?) v / cos~x. 



(4) 



This inequality gives the lower limit to the effective surface 
temperature as a function of %. It follows from Eq. (0) by 
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keeping the most important longitudinal thermal conduc- 
tivity but neglecting the transverse conductivity (which 
is often much lower). In this case the surface temperature 
T S (B, x) is naturally expressed through T» (B). If K|| » k± 
and re || cos 2 x^> K ± sin 2 X one can expect the real surface 
temperature to be very close to the lower limit. The in- 
equality (||) can be useful for testing numerical calcula- 
tions of T s in magnetized NSs. 

Equation (||) was first obtained by Greenstein & 
Hartke (1983). The authors derived (|J) as an estimate 
assuming actually that the thermal conductivities K|| and 
re_i_ are constant in the isolating layer. We emphasize that 
the equation is valid under much wider assumptions. Re- 
cently the solution similar to (||) has been proposed by 
Page (1995) for a dipole magnetic field. However, he used 
T S (B P ,0) instead of T\\(B), and T s (5 cq , tt/2) instead of 
T±(B) on the right-hand side (B p and B eq being the po- 
lar and equatorial magnetic fields, respectively). The ac- 
tual difference of his formula from Eq. (|J) is not large but 
Eq. (||) is thought to be more adequate. 

The advantage of Eq. (|[) is that it reduces the calcula- 
tion of T S (B, x) t° solving two heat conduction problems: 
for B normal (T\\(B)) and tangential (T±(B)) to the sur- 
face, as if at the magnetic pole and equator, respectively. 
The longitudinal (polar) problem has been considered very 
thoroughly by Van Riper (1989), and we use his results. 
We have taken the dependence of T s on T from his Fig. 29 
for B= 0, 10 11 , 10 115 , 10 12 , 10 12 - 5 , 10 13 , 10 135 G and for 
the surface gravity g s — 10 14 cm s -2 . We have interpo- 
lated these data for intermediate values of B. The trans- 
verse (equatorial) heat transfer problem has been analyzed 
by Schaaf (1991a). We have taken the dependence of T s on 
71 from his fitting Eqs. (29) and Table 3 for B= 10 u , 10 12 , 
and 10 13 G at g s — 10 14 cm s~ 2 . We have interpolated 
these data for intermediate values of B, and extrapolated 
them to lower B (if required) in such a way to reproduce 
the B = dependence of Van Riper (1989). We have used 
the scaling relation T s oc gl^ 4 to extend the results to the 
surface gravities of study (Sect. 3). 

Equation ([|) is valid for any magnetic field geometry. 
For simplicity, we assume that the surface magnetic field 
is dipolar. Let B p be the polar magnetic field, and let us 
introduce the reference frame which is flat for a local part 
of the NS subphotospheric layer. Taking into account the 
effects of General Relativity on the dipole field configura- 
tion (Ginzburg & Ozernoy 1964) one has 
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B R = B P cosO, B e = -B p f sin 9, (5) 
(1 - u) ln(l -u)+u - 0.5u 2 
* ~ ln(l - u) + m + 0.5m 2 ' 

where Br and Bq are, respectively, the radial and tangen- 
tial field components, is the colatitude (with O = at 
the pole), and u — R g / R. If u <C 1 the General Relativity 
effects are negligible and / = —1/2. With increasing u, 
the General Relativity effects make the field more radial. 
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Figure 1 shows the surface temperature distributions 
a NS with M = 1.3 Mq and R = 11.70 km for two 
■rnal temperatures: TJ = 10 s and 10 7 K, typical for 
neutrino and photon cooling stages, respectively, see 
ts. 3 and 4 below. The results can be explained as fol- 
s. The magnetic fields affect the heat transfer in two 



: M = 1.3 M Q R = 11.7 km : 

I T s = 10 s K : 




30 60 
Colatitudc (degrees) 



. 1. Distribution of the effective temperature T s from the 
(colatitude=0) to the equator (90) of the 1.3 Mq neutron 
with the dipole magnetic field for two values of the inter- 
temperature (Ti = 10 7 and 10 8 K) and several polar field 
ngths Bp. Open circles show the lower T s limit, Eq. (^), for 
= 10 12 G and T = 10 7 K. Thin long dashes — T a for the 
e dipolar field calculated without General Relativity effects 
= -1/2 in Eq. (5)). 

s (e.g., Yakovlev & Kaminker 1994). First, the fields 
> 10 10 G magnetize the electrons in the degenerate 
layer of the outer isolating layer (see above). Rapid 
tron Larmor rotation (ujbt > 1, where ojb is the gy- 
equency of electrons at the Fermi surface, and r is the 
tron relaxation time) suppresses the transverse con- 
tivity making the thermal conduction anisotropic. Sec- 
, the fields B > 10 12 G quantize electron motion - 
the structure of the electron Landau levels become pro- 
nounced in the thermodynamic and kinetic properties of 
degenerate and non-degenerate electrons. 

As seen from Fig. 1, the field B p < 10 10 G is too low to 
modify noticeably the surface temperature T s - the tem- 
perature is nearly isotropic over the surface. The magnetic 
fields 10 10 < B < 10 12 G magnetize but not quantize the 
electron motion. They strongly reduce the transverse ther- 
mal conductivity of degenerate electrons enhancing ther- 
mal isolation of the NS. These fields do not affect sig- 
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nificantly the longitudinal electron conductivity, the ra- 
diative thermal conductivity of outer non-degenerate sub- 
layers, and thermodynamic properties of matter. Accord- 
ingly the polar surface temperature Tji remains nearly the 
same as for B = but the equatorial temperature T± can 
be noticeably reduced. With increasing £?, the reduction 
becomes stronger, and T S (B, y) approaches the minimum 
value Tii (B)^/cos \ given by (Q) almost everywhere at the 
surface except very near to the equator. As a result, the 
total luminosity (0) decreases with growing B but this de- 
crease saturates after the inequality kh S> Kj_ is achieved. 

The quantizing magnetic fields B > 10 12 G induce 
strong quantum oscillations of the thermal conductivity 
of degenerate electrons. Generally, these fields increase the 
longitudinal electron conductivity. They also increase the 
longitudinal and transverse radiative thermal conductivity 
of non-degenerate layers. Moreover, the strongly quantiz- 
ing fields reduce the pressure of degenerate electrons. All 
these effects lower the thermal isolation. Accordingly, they 
increase the polar surface temperature while the equato- 
rial surface temperature remains strongly damped by the 
magnetization of degenerate electrons (Fig. 1, curves with 
B p = 10 13 G). The net effect is to enhance the total ther- 
mal luminosity L of the star. With increasing B, the effect 
does not saturate but becomes more pronounced. The ef- 
fects of the magnetic fields are stronger for a colder NS, 
where the isolating layer is thinner and more sensitive to 
the magnetic fields (e.g., Kaminker & Yakovlev 1994). 

Therefore the thermal luminosity is suppressed by the 
magnetizing fields in the equatorial region and enhanced 
by the quantizing fields in the polar region. Two effects, 
the suppression and the enhancement, act in different 
ranges of B and in opposite directions. The traditional 
calculations of NS cooling assume that the magnetic field 
is radial (polar) everywhere over the surface. Accordingly 
the calculations take into account the enhancement but 
neglect the suppression. The model surface temperature 
distribution (|J) allows one to include both effects at once. 

The validity of Eq. (|) can be checked by direct 2D 
simulations of the heat transport problem in the sur- 
face layers of a magnetized NS. The first attempt to per- 
form 2D simulations was made by Schaaf (1991b) under 
many simplified assumptions. Nevertheless, the 2D effects 
such as the heat flow from the hotter polar regions to 
the cooler equatorial ones produced by heat conduction 
or/and meridional and convective motions should only 
smooth the temperature variations over the NS surface. 
In this respect T s given by (||) can be considered as an 
upper limit of the possible temperature variation over the 
NS surface. It seems sufficient for our cooling calculations. 

3. Cooling calculations 

In order to analyze the effects of the magnetic fields on 
the NS cooling we have performed a series of numerical 
simulations. We have used the cooling code of Gnedin & 



Yakovlev (1993) and Gnedin et al. (1994) based on the 
approximation of isothermal NS interior. The cooling con- 
sists of two stages: the neutrino and photon ones. At the 
initial neutrino stage (i < (10 5 - 10 6 ) yrs) the NS cools 
mainly due to the emission of neutrinos from its interior. 
The internal temperature is ruled by the neutrino luminos- 
ity, and it is independent of the surface temperature; the 
surface thermal radiation just follows the internal cooling. 
At the later photon stage the star becomes colder, and the 
neutrino production mechanisms are inefficient. The star 
cools mainly via thermal surface radiation which governs 
the fall of the internal temperature. 

In our simulations, the equation of state of matter in 
the NS core has been taken from Prakash et al. (1988). Ac- 
cording to this equation of state, the core consists of neu- 
trons, protons and electrons (no hyperons or exotic par- 
ticles). The maximum mass of a stable NS is 1.7 M©. We 
have chosen two representative models, with M—1.3Mq 
and 1. 44 M . 

The model M=1.3M (i?=11.70 km, the central den- 
sity p c — 1.12 x 10 15 g cm~ 3 ) is typical for a NS with the 
standard neutrino energy losses. In the cores of such stars, 
neutrinos are produced by the modified Urea processes 
and the nucleon-nucleon neutrino-pair bremsstrahlung 
(Friman & Maxwell 1979, Yakovlev & Levenfish 1995). 
The modified Urea reactions are n + N — ► p + e + N + D e , 
p+e+N — ► n+N+v e , where A is a spectacular nucleon re- 
quired to satisfy momentum conservation. Until recently, 
it has been widely assumed that the main contribution 
into the modified Urea comes from the neutron branch 
reaction (A~ = n). However according to Yakovlev & Lev- 
enfish (1995) the proton branch reaction (A — p) is of 
comparable efficiency, and we have included it into the cal- 
culation. The nucleon-nucleon bremsstrahlung (A + A" — > 
N + N + v + v) consists of three branches (n + n, n + p, 
p + p), and all of them have been incorporated using the 
results of Friman & Maxwell (1979) and Yakovlev & Lev- 
enfish (1995). The effective masses of neutrons and pro- 
tons (which enter the expressions for the neutrino energy 
loss rates and heat capacities) have been set equal to 0.7 
of the bare masses. 

The model M=1.44M (#=11.37 km, p c = 1.37xl0 15 
g cm -3 ) represents a massive NS whose neutrino luminos- 
ity is greatly enhanced (Lattimer et al. 1991) by the direct 
Urea process (n — > p + e + D e , p + e — > n + v e ). The latter 
process intensifies the neutrino energy production rate by 
about 5 orders of magnitude for T\ ~ 10 9 K. The direct 
Urea is a threshold reaction: it is allowed (Lattimer et al. 
1991) for certain equations of state at rather high densities 
(to satisfy momentum conservation of the reacting parti- 
cles). For the given equation of state, the critical threshold 
density is p cr = 1.30 x 10 15 g cm~ 3 . In our NS model, the 
direct Urea is allowed in a small central kernel of the star 
with p > p c ; the kernel radius is 2.32 km and its mass is 
0.035 Mq. The neutrino luminosity of the kernel greatly 
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exceeds the neutrino luminosity from all other parts of the 
star. 

In both NS models we have also taken into account 
the neutrino luminosity of the stellar crust produced by 
neutrino-pair bremsstrahlung of electrons colliding with 
atomic nuclei. This has been done in the same simpli- 
fied manner as in the work of Maxwell (1979). We have 
neglected the contribution of the crust into the NS heat 
capacity, since the crust mass is very low in our models. 
The mass of the crust in a NS with stiff equation of state 
and/or in a low-mass NS can be larger, i.e., in principle, 
the crust can contribute noticeably into the heat capacity. 
In our models, the heat capacity is produced by neutrons, 
protons and electrons in the NS core. We neglect possible 
superfluidity of neutrons and protons since our main goal 



is to study the influence of the surface magnetic fields on 
NS cooling. The latter influence is determined by isolating 
properties of the NS surface layers (Sect. 2), and it is not 
connected directly with the internal superfluidity. 

We have adopted the dipole magnetic field geometry, 
Eq. (5), as described in Sec. 2. Using Eqs. (||) and (5) 
we can easily calculate the photon thermal NS luminos- 
ity which is an important ingredient of the cooling code. 
We have checked that the General Relativity effects in 
the magnetic field distribution (5) are almost negligible 
(Fig. 1) for our NS models. For comparison, we have also 
considered the case when the magnetic field is radial and 
has the same strength everywhere over the NS surface. 
This monopolar field geometry is traditional for the NS 
cooling theories. The field strengths have been varied from 
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B = to B ~ 10 13 G. We have followed the cooling until 
the mean effective surface temperature drops below 10 5 K. 
It has been unreasonable to proceed further since the re- 
sults of Van Riper (1989) and Schaaf (1991a) used for 
calculating T s are obtained for a not too cold NS. 

4. Results and discussion 

Figures 2 and 3 show the standard (1.3 Mq) cooling curves^ 
for the dipolar and monopolar magnetic fields of different^ 
strengths. The curves display redshifted NS photon lumi- 
nosity versus age t. Since the surface temperature is dis- 
tributed nonuniformly, the NS radiation, as detected by a 
distant observer, depends on angle between line of sight 
and the magnetic axis. We present the luminosity from 
the total NS surface which is important for NS cooling 
theories but does not necessarily determine the radiation 
flux in a specific direction. The field geometry is seen to 
be very important. The monopolar magnetic field always 
reduces the NS thermal isolation (only the electron quan- 
tization effects operate). Accordingly, the monopolar field 
increases the photon luminosity L at the neutrino cool- 
ing stage (t < 3 x 10 5 ) K, and decreases L at the subse- 
quent photon cooling stage. This is typical for the tradi- 
tional cooling calculations (e.g., Van Riper 1991). On the 
contrary, the dipolar field B p = 10 11 G, for example, de- 
creases L at the neutrino cooling stage, and increases at 
the photon stage. This is certainly because the magnetiz- 
ing fields enhance the equatorial thermal isolation. How- 
ever, for larger B p , the polar drop of the thermal isolation 
becomes more important and produces just the traditional 
cooling effect: it raises L at the neutrino cooling stage and 
decreases at the photon stage. With increasing B p , the po- 
lar effect becomes stronger than the equatorial one, but 
both effects interfere, and the results are quantitatively 
different from the traditional ones. 

The above effects are also illustrated in Fig. 4 where 
we plot L(B)/L(0) versus B for the 1.3 M Q NS of sev- 
eral ages. For the monopolar field, the dependence of 
L(B)/L(0) on B at given age is monotonic, and one needs 
B > 10 12 G to affect the cooling. For the dipolar field, 
the same dependence is non-monotonic, and the fields 
B > 10 10 G affect noticeably the cooling curves. It is re- 
markable that the dipolar fields B p « 3 x 10 12 almost 
do not affect the NS cooling. For these fields, the en- 
hancement of the equatorial isolation nearly compensates 
the drop of the polar one, and the cooling proceeds just 
as if the NS were non-magnetized. Nevertheless the field 
B p m 3 x 10 12 makes the surface temperature distribu- 
tion very anisotropic (Fig. 1), which can produce strong 
modulation of the surface thermal radiation (Shibanov et 
al. 1995). When B p <3x 10 12 G, the equatorial effects 
prevail while for B p > 3 x 10 12 G the traditional po- 
lar effects are stronger. Note that strong dependence of 
L(B) / L(0) on B for lg(t[yr]) = 6 and 6.3 is partially ex- 
plained by steep slopes of the cooling curves (Fig. 2) at 
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Fig. 4. Ratio L(B) / L(0) of the photon luminosities of magne- 
tized and non-magnetized 1.3 Mq neutron stars with dipolar 
(Fig. 2) and monopolar (Fig. 3) fields versus field strengths for 
several stellar ages t[yrs] (values of lg(t) are given near the 
curves) 



the photon cooling stage. Figures 5 and 6 plot the cooling 
curves for the 1.44 Mq NS. The curves for this (rapidly 
cooling) NS are qualitatively the same as for the standard 
cooling (Figs. 2 and 3). However the photon luminosity L 
is much lower than for the standard cooling since the rapid 
cooling goes much faster at the neutrino cooling stage. The 
appropriate luminosity ratios L(B) / L(0) versus B for sev- 
eral stellar ages are shown in Fig. 7. The curves are quite 
similar to those for the standard cooling (Fig. 4). 



5. Conclusions 

We have used a realistic model (||) of the anisotropic dis- 
tribution of the effective temperature T s over the surface 
of a magnetized NS and have calculated NS cooling. We 
have considered both standard and rapid cooling of the NS 
with the dipole magnetic field. The results are noticeably 
different from the traditional results obtained for radial 
(monopolar) magnetic fields. We have shown (Figs. 2-7) 
that the dipole fields 10 10 <B p <3x 10 12 G decrease the 
stellar photon luminosity at the neutrino cooling stage (by 
a factor of ~ 2), and increase the luminosity at the pho- 
ton cooling stage (by about one order of magnitude), as 
compared to a non-magnetized NS, contrary to the tradi- 
tional theories developed for monopolar magnetic fields. 
The effect is produced by the growth of the thermal isola- 
tion of the equatorial surface layers due to the low thermal 
conductivity of degenerate electrons across the magnetic 
field. On the other hand, the dipole fields B p > 3 x 10 12 G 
increase the photon thermal luminosity at the neutrino 
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Fig. 6. Same as in Fig. 3 for the 1.44 Mq neutron star. 



cooling stage and decrease the luminosity at the photon 
cooling stage, in comparison with the luminosity at B = 0. 
This conclusion agrees qualitatively with the results of the 
traditional theories but the actual magnetic field effect is 
much weaker than predicted by these theories due to the 
interference with the opposite effect in the equatorial re- 
gions. We have obtained also that the dipole field with 
B p »3x 10 12 G has almost no effect on the NS cooling. 

Our results are based on a model of the surface temper- 
ature distribution applied to the dipole magnetic field. It 
would be easy to consider the NS cooling for other mag- 
netic field geometries (quadrupolc, shifted dipole, etc.), 
and the results are expected to be qualitatively the same 
as for the dipole field: the parts of the surface with essen- 



tially non-radial magnetic fields should noticeably increase 
the thermal isolation (Sect. 2). 

For further studies of the magnetic field effects on the 
NS cooling, it would be highly desirable to perform de- 
tailed investigation of the heat transport in the outer NS 
layers using the best microscopic physics available (equa- 
tion of state, thermal conductivities, etc), and obtain thus 
exact surface temperature distribution in magnetized NSs. 
The distribution (^|) used above is expected to give an up- 
per limit for the surface temperature variation produced 
by the magnetic fields (Sect. 2). It can be used for a test 
in more advanced theories. Even if (|§|) is not very accu- 
rate at large magnetic fields near the equatorial regions, 
it can yield quite accurate results since the cold equato- 
rial surface parts do not contribute significantly into the 
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Fig. 7. Same as in Fig. 4 for the 1.44 Mq neutron star. 

useful for theoretical interpretation of thermal X-ray ra- 
diation of NSs. All the objects from which this radiation 
has been detected (Sect. 1) are thought to possess mag- 
netic fields of about several times of 10 12 G. We can expect 
that such fields do not affect significantly the NS cooling 
but they can produce strongly anisotropic surface temper- 
ature variation and associated modulation of the surface 
thermal radiation. 

The magnetic fields broaden the allowed ranges of 
the photon luminosities in the L-t diagram for NSs of 
various masses, radii and equations of state with stan- 
dard or enhanced neutrino energy losses (e.g., Ogelman 
1994). According to the traditional theories, high mag- 
netic fields increase the upper boundary of the expected 
photon luminosity at the neutrino cooling stage (typically, 
at 10 2 < t < 3 x 10 5 yrs) and decrease the lower bound- 
ary of L at the photon stage. Our results (Figs. 4 and 
7) show that the dipolar fields B ~ 10 11 G decrease the 
lower boundary of L at the neutrino stage and increase 
the upper boundary at the photon stage. 
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